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ABSTRACT 

The centre of our Galaxy harbours a 4 mihion solar mass black hole that is unusually quiet: its 
present X-ray luminosity is more than 10 orders of magnitude less than its Eddington luminosity. 
The observation of iron fluorescence and hard X-ray emission from some of the massive molecular 
clouds surrounding the Galactic Centre has been interpreted as an echo of a past 10'^^ erg s~^ flare. 
Alternatively, low-energy cosmic rays propagating inside the clouds might account for the observed 
emission, through inverse bremsstrahlung of low energy ions or bremsstrahlung emission of low energy 
electrons. Here we report the observation of a clear decay of the hard X-ray emission from the 
molecular cloud Sgr B2 during the past 7 years thanks to more than 20 Ms of INTEGRAL exposure. 
This confirms the decay previously observed comparing the 6.4 keV line fluxes measured by various 
X-ray instruments, but without intercalibration effects. The measured decay time is 8.2 ± 1.7 years, 
compatible with the light crossing time of the molecular cloud core . Such a short timescale rules out 
inverse bremsstrahlung by cosmic-ray ions as the origin of the X ray emission. We also obtained 2-100 
keV broadband X-ray spectra by combining INTEGRAL and XMM-Newton data and compared them 
with detailed models of X-ray emission due to irradiation of molecular gas by (i) low-energy cosmic- 
ray electrons and (ii) hard X-rays. Both models can reproduce the data equally well, but the time 
variability constraints and the huge cosmic ray electron luminosity required to explain the observed 
hard X-ray emission strongly favor the scenario in which the diffuse emission of Sgr B2 is scattered 
and reprocessed radiation emitted in the past by Sgr A*. The spectral index of the illuminating 
power-law source is found to be F '^ 2 and its luminosity 1.5 — 5 x 10^^ erg s~^, depending on the 
relative positions of Sgr B2 and Sgr A*. Using recent parallax measurements that place Sgr B2 in 
front of Sgr A*, we find that the period of intense activity of Sgr A* ended between 75 and 155 years 
ago. 
Subject headings: Galaxy : center, ISM: clouds (Sgr B2), X-ray : ISM 



1. INTRODUCTION 

The supermassive black hole o f 4.3 x 10^ Mq 
(jGhez et al.l 120081: iGillessen et"all I2OO90 that sits at the 
centre of the Galaxy is remarkably quiescent nowadays 
with an X-ray luminosity of about 10^'^"'^''^ ergs~^, more 
than 10 orders of magnitude less than its Eddington 
luminosity. The brightest X-ray flare rec orded so far 
reach ed a level of a few 10'^^ ergs~^ (jPorauet et al.l 
l2008f ). Active Galactic N uclei have a relatively low duty 
cycle (-- 10"^ . [HoI 120081 ). so that they spend most of the 
time in a low luminosity state. It is therefore not ex- 
cluded that Sgr A* has experienced periods of increased 
activity in the past. Looking for evidence of such a past 
activity of Sgr A*, we turn to its surroundings that 
would probably carry the scars of suc h an event. 

Th e central molecular zone or CMZ ([Morris fc Serabvnl 
|1996[ ). although relatively small in size with a radius of 
about 200 pc, contains about 10% of the Galaxy's molec- 
ular mass. It is populated by some of the most massive 
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(10'' - lO*" Mq) and dense (> 10^ cm"^) giant molecu- 
lar complexes in the Galaxy. Sgr B2 is the densest and 
most massive among them with a mass of a few IO^Mq, 
and density in its 5 pc core reaching 10^ cm~"^. It has an 
envelope that e xtends to ~20 pc with a n average density 
of - 10^ cm-3 (jLis fc Goldsmith|[l99l . 

The CMZ is also a source of non-thermal radiation. 
ART-P on board GRANAT was t he first to detect har d 
X-ray emission from the region fSunva ev et al.l I1993[) . 
The ASCA satellite discovered a strong iron Ka fluo- 
rescence line at 6.4 k eV with a large equivalent width of 
1-2 keV from Sgr B2 (jKovama et al.|[l996i) . Chandra has 
surveyed the central region of our Galaxy and mapped 
the dis tribution of Fe li ne emission from the molecular 
clouds (jPark et al.ll2004[ ). Today, thanks to the improved 
sensitivity of INTEGRAL/IBIS above 20 keV, the emis- 
sion of Sgr B2 was res olved from the rest of t he Galactic 
Centre (GC) sources (jRevnivtsev et al.|[2003 ). 

The features seen in Sgr B2 are not unique 
among clouds in the CMZ. GO. 13-0. 13 (also konwn 
as GO. 11-0. 11) displays strong 6.4 keV line emission 
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(|Yusef-Zadeh et al]|2002[ ). as well as a hard X-ray c ontin- 



uum detected at energies above 50 keV (Bclangcr et al 
2006 ) . Sgr C shows Fe Ka line emission (jNakaiinia et al 



20091 ) ■ however its associated hard X-ray continuum 
cannot be resolved by current instruments because 
it lies close to a br ight X-ray binary, KS 1741-293 
(|de Cesare et al.ll2Q07l ). 

The origin of this non-thermal emission is contro- 
versial. A likely explanation for the continuum and 
the strong Ka line, is Compton scattering and K- 
shell photo-ionization of neutral or low-ionized iron 
atoms by ai i intense X-ray rad i ation field generated 
in the GC (iSunvaey et all Il993t iKovama et all 119961: 
iSunvaev k, Churazovlll998[ ). Here, the clouds would act 
as X-Ray Reflection Nebulae (XRN). Since there has 
been no report of a sufficiently bright X-ray transient 
(L > 10"^^ erg/s) close to Sgr B2 in the last 20 years, it 
has been argued that the origin of the emission should 
lie in the most central regions and have a luminosity of 
at least 10^^ erg/s for several years. The most probable 
candidate is Sgr A*. This, however, requires an increase 
of more than 6 orders of magnitude compared to its cur- 
rent luminosity (Koyama ct al. 1996). It has also been 
suggested that the Sgr A East supernova remnant in- 
teracting with the GC massive molecular clouds could 
provide the required luminositv ()Frver et al.|[2006l ). 

The propagation of energetic charged particles inside 
the molecular clouds has also been invoked as an explana- 
tion of the observed X-rays from the region. For instance, 
lYusef-Zadeh et al.l 12007.) suggested that the density of 
cosmic rays (CR) is larger in the CMZ than in the rest of 
the Galaxy. Their hypothesis is based on several results: 
the observation of diff use low-frequency r adio emission 
in the central regions ( LaRosa et al.l 120051 ). the high es- 
timates of the ionization rat e compared to t he values 
obtained in the Galactic disk (|Oka et al.ll2005l ). and the 
detection of hard interst ellar emission from the molecular 
clouds at TeV energies (jAharonian et al.ll2006[ ). 

Two distinct processes can produce hard X-ray contin- 
uum and neutral iron line emissions. First, low-energy 
cosmic-ray electrons (LECRE), with energies up to a few 
hundred keV, diffusing into dense neutral matter will 
produce nonther mal bremsstra.hlung and collisional K- 
shell ionization (jValinia et al.l I2OOO0 . The presence of 
strong non-thermal radio filaments in the vicinity of sev- 
eral molecular clouds makes it reasonable to consider the 
presence of a large flux of c osmic ray electrons in th ese 
regions see e.g. GO.13-0.13 (lYusef-Zadeh et al.l l2002|) or 
Sgr C (jYusef-Zadeh et al.ll2007D . Second, subrelativis- 
tic ions propagating in the molecular cloud can radi- 
ate through inverse bremsstrahlung in the hard X-ray 
domain and create iron Ka vacancies. Knock-on elec- 
trons generated by collisions of primary ions with the 
gas also cont ribute significantly to the total emission 
(jDogiel et al.l [20091. These subrelativistic ions can also 
be invoked to explain the putative hot plasma pervading 
t he GC. 

iRevnivtsev et all (2004) have argued that the X-ray 
emission from Sgr B2 is difficult to explain by LECREs, 
because the large power that cosmic-ray particles would 
need to deposit in the molecular gas to account for 
the X-ray flux would be comparable to the bolometric 
(mostly infrared) luminosity of the complex. They also 
put into question the large metallicities implied by the 



large equivalent width of the line produced by collisional 
K-shell interaction. 

To help resolve this ambiguity, Sgr B2 and the other 
clouds of the CMZ have been the subject of numerous 
observation campaigns. In particular, a deep Chandra 
observation revealed tha t the 6.4 keV line em ission is 
oriented towards the GC (jMurakami et al.ll2001[ ). It also 
allowed the contribution from individual sources in the 
CMC associated with dust condensations (e.g. Sgr B2 
North (N) Middle (M) and South(S) to be distinguished 
^Takagi et al.J 2002). 

Recently, thanks to the monitoring of the region by 
many instrume nts, time dep e ndent effects have been ob- 
served. First, iMuno et al.l (|2007[ ) have detected time 
varying diffuse emission in the 4-8 keV range a few ar- 
cmin away from Sgr A* comparing Chandra images taken 
before and after 2004. Since this region is emitting the 
6.4 keV line, it is likely that the varying diffuse emission 
can be linked with iron fluorescence, making this cloud 
a likely X-ray reflection nebula. More recently, com- 
paring recent observations with various X-ray telescopes 
(namely Suzaku, XMM-Newton, Chandra), llnui et aD 
(2009) were able to show a temporal decrease in the Ka 
line emission from Sgr B2, strengthening the reflection 
interpretation for this cloud as well. 

The hard X-ray continuum carries many clues as re- 
gards the origin of the emission: its spectral shape and 
its temporal behaviour in particular. With its sensitivity 
and with its very long exposures of the Galactic central 
region, INTEGRAL can therefore provide invaluable in- 
sight into the nature of the hard X-ray emission of the 
CMZ. We present here the lightcurve and spectrum of the 
Sgr B2 emission taken over 7 years of observations and 
show that the Sgr B2 flux has decreased by more than 
40% since 2003. We then study the combined XMM- 
Newton/INTEGRAL spectrum obtained from observa- 
tions prior to 2004 which we use to test LECRE and 
XRN models. Finally, we discuss the implications of our 
results in relation to the likely nature of the non thermal 
emission from the GC molecular clouds. 

2. OBSERVATIONS AND DATA ANALYSIS 
2.1. IBIS/ISGRI data 

Since its launch in October 2002, INTEGRAL 
(jWinkler et al.l '2003) has been performing deep obser- 
vations of the Galactic Centre region for more than 1 Ms 
every 6 months, first as part of the core program Galactic 
Centre Deep Exposure and then with the Galactic Cen- 
tre Key Programme. A number of dedicated pointings 
as well as Target of Opportunity observations also add 
exposure to the GC thanks to the wide field of view of 
the main instruments on board INTEGRAL. We have 
therefore used all public data available within 13° of the 
GC source position, as well as the latest observations 
of the Galactic Centre Key Progamme. We applied a 
number of selection criteria on the science windows' (i.e. 
individual pointings) quality so as to improve the overall 
dataset quality (imposing an exposure larger than 600s 
per pointing and rejecting images with bad significance 
distribution). This yields a number of 7600 science win- 
dows for a total of about 20 Ms of exposure. 

The size of the point spread function (PSF) is of funda- 
mental importance when determining the contribution of 



Fading hard X-ray emission from Sgr B2 



1.0001 



0.500 





























IGRJ1 7456-2901 



T5 
3 


0.000 


XI 




+- 




a 




o 








u 




ra 






-0.500 



IGR J1 7475-2822 ^^ 1743^284^^ ^^ 1741-29? 



IE 1740.7-2942 



IGR J1 7497-2821 SAX J1 747.0-2853 



1 A 1742-294 



-1.000 






LX 1744-299 



1.500 1.000 0.500 0.000 359.500 

Galactic longitude 



359.000 358.500 



16:00.0 



Fig. 1. — IBIS/ISGRI significance map of the Galactic central regions in the 20-40 keV energy range. The image was produced combining 
all the GC observations. The main sources are indicated by the black crosses. The position of IGR J17475-2822 is compatible with that of 
the Sgr B2 molecular cloud. 

ages of all science windows in 3 energy bands: 20-40 
keV, 40-60 keV and 60-100 keV. Individual images were 
then mosaicked. A close-up of the Galactic central de- 
grees is visible in Fig [TJ In order to extract regularly 
sampled lightcurves, we divided the observations into 14 
periods (2 per year), each corresponding to an observa- 
tion window lasting up to one month around September 
and March of each year, and we produced the 14 images 
of the Galactic Centre spanning from March 2003 up to 
September 2009. 

Imperfect background and detector uniformity correc- 
tions introduce systematics biases and errors in the fi- 
nal mosaic images. The large number of bright sources 
in the Galactic central region also creates important 
systematic effects because of the remaining uncertain- 
ties in the IBIS/ISGRI PSF model that is subtracted 
from each individual h xiage to remove ghost sidelobes 
(jGoldwurm et al.ll2003[) . To take into account all these 
systematics in the final mosaics, we applied a correction 
factor to the variance skymaps to recover a normalized 
significance distribution of each image. 

To extract count rates of sources in the central degree 
where source confusion is large, we have to fit simulta- 
neously the fiuxes of all of the 11 sources present there. 
Since Sgr B2 is in a somewhat less crowded region, it is 
sufficient in this case to model the fiuxes of just three 
point sources (namely IE 1743.1-2843, IGR J17497-2821 
and IGR J17475-2822/Sgr B2) using a gaussian PSF of 
width a— 0.1° in a 1° wide square region around the po- 
sition of Sgr B2. 

Since there is a correlation from bin to bin, due to the 
correlation technique used to make images with coded 
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Fig. 2.— XMM-Newton/MOS 6.4 keV line flux image of the Sgr 
B2 field. The cyan circle gives the best fit position of IGR J17475- 
2822 with the associated 90% confidence error on the position. The 
white dashed circle shows the 4.5' wide extraction region used for 
the combined spectral analysis. 



the various sources in the GC. We therefor e focused only 
on th e IBIS/ISGRI data with its 12' PSF (jLebrun et al.l 
l2003f ). We analyzed the whole dataset using the Off- 
line Science Analysis (OSA) software version 7.0 with 
an improved energy reconstruction taking into account 
gain losses during the whole mission. We produced im- 
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Fig. 3. — The central degree as seen by IBIS/ISGRI on board INTEGRAL in the 20 to 60 keV energy range. Each row corresponds to a 
year from 2003 to 2009. The many sources in the GC display strong time dependence. The green circle shown in the top panel corresponds 
to the position of the core region of the Sgr B2 cloud -see the text. The associated hard X-ray source, IGR J17475-2822, shows a clear 
decline in flux after the 7 years of monitoring. 



mask instruments, we cannot use the errors computed 
with the chi-square fit for the source fluxes. We therefore 
used the corrected image variance as the error on the 
source fluxes, which is a conservative estimate. 

The ISGRI camera suffers from efficiency losses with 
time. To take these into account, we performed the same 
analysis on the calibration observations of the Crab Neb- 
ula taken on-axis during each of the 14 observation win- 
dows on the GC. We got a lightcurve of the measured 
Crab count rate in each energy band showing a clear de- 
crease (about 5% during the whole mission in the 20-40 
keV energy range), which is well fitted by a linear de- 



cay. This provided us a time-dependent correction to 
the measured count rates which we applied to the GC 
sources. As a verification procedure, we used the Ophi- 
uchus cluster as a secondary calibrator. It is visible in all 
GC mosaics and is exte nded with respect to the ISGRI 
PSF (jEckert et al.ll2008D . We extracted its flux in each of 
the 14 mosaics and found the resulting lightcurve, after 
application of the flux correction, to be compatible with 
a constant. 



2.2. XMM-Newton data 
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Sgr B2 was observed by XMM-Newton for 50 ks in 
September 2004 (Obsid: 0203930101). Event files were 
cleaned for flaring events and background-subtracted 
mosaicked images of the EPIC/MOS instruments have 
been produced. The background used for the subtrac- 
tion is taken from reference high-latitude observations 
obtained in the same observing mode and filter and 
includes instrumenta l as well as cosmic backgrounds 
(jCarter fc Readl2007f ). Prior to its subtraction, the back- 
ground was normalized using the measured flux in the 10- 
12 keV energy range. To correct for the exposure and the 
vignetting, the mosaicked background-subtracted count 
images were divided by their mosaicked exposure maps 
which include the vignetting corrections. The adaptative 
smoothing tool asmooth was applied to the images. 

Images were produced in the 4-6 keV and in the 6.3 - 
6.48 keV energy ranges. To remove the continuum from 
the 6.4 keV line image, we used the background sub- 
tracted 4-6 keV image, renormalized according to the 
fitted shape of the overall spectrum in the field of view. 

We extracted spectra only from the 2004 observation, 
which was taken after the INTEGRAL launch and can 
be used for combined spectral fitting. We used a 4.5' 
radius extraction region. Because of the large spectral 
extraction region radius, we extracted background events 
from empty field observati ons taken in a similar time 
frame (|Carter fc Readl[2007[ ). 

3. RESULTS 
3.1. Time-varying hard X-ray emission from Sgr B2 

In the complete mosaic (covering 7 years of data) , the 
source IGR J17475-2822, associated with Sgr B2, is de- 
tected at a confidence level of more than 50 cr; see Fig 
[TJ We left its position free in the multi-source fitting 
procedure described in section 2.1 and found it to be 
1 = 0.67°, b = -0.05° with a 90% confidence radius of 
0.7' (|Gros et al.l I2003D . in perfect agreement with the 
Sgr B2 core (MO. 66-0. 02) position, see the comparison 
with the 2004 XMM-Newton image in the Fe Ka band 
(figure [5]). Nevertheless, the large IBIS PSF encloses a 
significant part of the Sgr B2 complex, so that a contri- 
bution from neighbouring structures such as MO. 74-0. 09 
cannot be excluded. In particular, the larger and more 
diffuse component of the molecular cloud should be con- 
tributing to the hard X-ray fiux. But its flux should be 
diluted because of the extent relative to the PSF size. 
If we ass ume that the envelope diameter is of the order 
of 45 pc (|Lis fc Goldsmithlll990D . the expected diameter 
at 8 kpc is 20', significantly larger than the IBIS PSF. 
This implies a reconstructed peak fl ux of about 40 % of 
its equivalent point source flux (Ren aud et al.l[2006l ). 

Using the 2004 XMM-Newton/EPlC observation of the 
region, we checked that no point source in the field is 
bright enough to produce such a strong emission, so that 
we can unambiguously associate the hard X-ray source 
with diffuse emission from the Sgr B2 cloud. We also 
fitted the INTEGRAL source position at different time 
periods during the 7 years of observation but found no 
significant displacement of the source position. 

Looking at INTEGRAL images of the GC we produced 
on a year timescale, most of the sources in the Galactic 
central degrees show strong variability on shorter time 
scales, as was noted bv .Kuulkcrs et al.. (,2007,); see Fig 



[31 One can also note a clear but slow decrease of the 
fiux of the source compatible with Sgr B2. To properly 
quantify this effect, we built a lightcurve of the Sgr B2 
flux over the 7 years of INTEGRAL observations. To 
do so, we flxed the IGR J17475-2822 source position to 
the one obtained in section 3.1 and extracted its count- 
rate in each of the 14 time periods using the multi-source 
fitting procedure. To correct for the varying instrument 
response, we applied the count-rate correction described 
in section 2.1. The resulting lightcurve is shown on figure 
m It is the most regularly sampled lightcurve of the non- 
thermal emission from this molecular cloud obtained so 
far, with a data point every six months. It is also the 
only one in the hard X-ray domain. We com pare it to 
the 6.4 keV measurements of llnui et al.l (|2009[ ). 

The Sgr B2 20-60 keV emission shows indications of 
a count rate decrease since 2003, reaching 40% in 2009 
compared to the 2003 fiux. To estimate the significance 
of the effect, we performed an F-test comparing the x^ 
of a constant fiux fit (x^ = 40.06 d.o.f. = 13) and a 
linearly varying one {x^ = 15.16 d.o.f. — 12). The vary- 
ing flux is preferred at the 4.80" level over a flat light 
curve. The decay time (time to divide the flux by two) 
is found to be T = (4.0 ± 0.8) x 10^ days, that is, 8.2 
± 1.7 yrs, which is compatible with the light crossing 
time of the molecular cloud core and is compatibl e with 
previous expectations f Sunvaev fc Churazovl ll998'). The 
a mplitude of t he va riation is consistent with that found 
bv llnui et all (pOOl . 

3.2. Broadband spectral analysis 

Sgr B2 is a faint object located in a crowded region 
containing more than a hundred sources within the IBIS 
field. Therefore its spectrum cannot be extracted using 
standard OSA spectral routines. We have built mosaics 
in 15 energy bands ranging from 18 to 300 keV, grouping 
the data in 2 year intervals. We used a fitting procedure 
similar to that described in section 2.1 to fit the flux of 
Sgr B2 as well as that of neighboring sources and the 
background in each energy band. The position of Sgr B2 
was fixed to the position fitted in the complete dataset 
in the 20-60 keV energy range. We applied the same 
procedure to estimate the errors on the fiux measurement 
to take into account systematic errors due to imperfect 
background and source side lobe subtraction. 

Because of the evolution of the instrument response 
with time, we computed mean ancillary response func- 
tions, arf^ for each of the 2-year periods using matrices 
distributed with the OSA. Applying these to the Crab 
Nebula observations, we estimated a systematic uncer- 
tainty of ^ 8%, which we applied to all spectral fits of 
the IBIS/ISGRI data. 

We performed a broadband fit to the 2004 XMM- 
Newton EPIC/MOS and PN and the 2003-2004 
IBIS/ISG RI spectra. We f ollowed an approach similar 
to that of iRvu et al.l ()2009[ ) but considered only events 
above 3.5 keV to avoid modeling the low energy back- 
ground in the GC and did not take into account the 
Cosmic X-ray Background (CXB) contribution which is 
subtracted by the blank sky events. We therefore used 
an APEC model for the hot thermal emission (the so- 
called 8 keV plasma) and we modeled the spectrum from 
the molecular cloud as the sum of a power law and of 
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Fig. 4. — (Black circles) Light curve of the corrected 20-60 keV flux in mcrab of Sgr B2 as measured with IBIS/ISGRI from 2003 to 
2009. The hght curve of the secondary calibrator (Ophiuchus cluster) obtaine d in the same con ditions is shown for comparison {Light blue 
squares). Superimposed {red squares) are the Fe Ka line fluxes obtained bv llnui et al.l |2003) with ASCA, Chandra, XMM-Newton and 
Suzaku observations taken from 1994 to 2006. The hard X-ray flux shows a clear decay of up to 40% of the measured Sgr B2 flux from 
2003 to 2009. A chi-square hypothesis test favors a linearly decreasing flux hypothesis over a constant flux at a 4. Sit level, the line shows 
the best flt. The 6.4 keV measurements show a similar trend after the year 2000, the flux being apparently constant from 1994 to 2000. 



two gaussians for the iron Ka and K/3 lines (respec- 
tively at 6.4 and 7.06 keV). All the emission components 
were subject to line of sight absorption due to a fore- 
ground interst ellar column density of ~ 6 x 10^^ cm~^ 
(|Sakanoet alll2002| ). The resulting XSPEC model is 
therefore wabs x (apec -I- gaus + gaus -I- wabs x pegpw) . 
The energies of the two iron lines are fixed to 6.4 and 
7.06 keV respectively, and we impose the ratio of F e K/3 
to Ka to be equal to 15% (iMurak ami et al.ll2000D . We 
obtain a satisfactory fit with a reduced x^ of 1.19. The 
plasma temperature is 5.8 ± 0.3 keV and the power-law 
component has a spectral index of 2.07 ± 0.05 and a 
2-100 keV unabsorbed flux normalization of 8.4 ± 0.4 
xlO^^^erg cm~^ s~^. The flux of the iron Ka line is 
found to be 8.9 ± 0.6 xlO^^ cm^^ s~^ . The intrinsic 
absorption is found to be 6.8 ± 0.5 xlO^^ cm~^, which is 
very large but still lower tha n estimates of the maximum 
column density of the cloud (jProtheroe et al.| [2008). The 
foreground absorption is 6 .9 zb 0.5 xlO^^ cm~^, consis- 
tent with other estimates ()Rvu et al.ll2009D . 

In the following, we try to apply detailed emission mod- 
els to the Sgr B2 spectrum to constrain the origin of the 
emission. 

4. ON THE NATURE OF THE HARD X-RAY EMISSION 

Our observation of a clear decay of the hard X-ray flux 
is consistent wit h that reported earlier for the 6.4 keV 
line emission by llnui et all ()2009[ ) - see Fig 01 We mea- 
sure a timescale of 8.2 ± 1.7 yrs for a 40% linear decrease 
of the hard flux compared to what was measured in 2003. 



The observation of such a rapidly varying fl ux is very dif- 
ficult to explain with cosmic ray particles. iDogiel et al.l 
(2009) have invoked a population of 100 MeV protons as 
a possible origin of the hard X-ray emission from the GC 
clouds. But since the relevant Coulomb cooling time, 
even in a dense cloud such as Sgr B2, is of the order 
of several hundred years, we can exclude the possibility 
that such a component makes a significant contribution 
to the observed X-ray emission. 

High energy {^ GeV) electrons can also produce hard 
X-ray emission through inverse Compton scattering on 
the ambient IR emission which is expected to be high 
in the Sgr B2 c ore b ecause of its huge luminosity. 
iGoicoechea et all ()2004D found a far IR luminosity of 
8.5 X 10^ L0 originating mainly from the Sgr B2(M) and 
(N) cores. In a radius of 1 pc, this creates a huge 
radiation density of 6 x lO^eV cm~^ in FIR photons. 
iLeahv fc Creightoiil (|1993l ) show that for a 1/r^ density 
profile, the Fe Ka EW can reach about 1 keV for an 
optical depth of the order of 1. With column densi- 
ties towards the core of the cloud reaching 10^** cm~^ 
(jProtheroe et al.l 120081 ). inverse Compton emission of 
high energy electrons in the centre of Sgr B2 might cre- 
ate a signal similar to that which we observe. Neverthe- 
less we are faced with similar problems: with densities 
reaching 10^ cm"'^ the dominant energy loss mechanism 
is bremsstrahlung, but its characteristic timescale is of 
the order of 50 yrs. 

We are left only with low energy cosmic ray electron 
models. In this case, variability of the overall flux is not 
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excluded since the energy loss timescale for ionization 
losses of ~ 100 keV electrons is of the order of or less than 
1 yr in dense environments of 10^ cm~^. We discuss this 
possibility in more detail in the next subsection using a 
detailed spectral model of the emission. 

4.1. A spectral model of LECRE emission 

In order to test for a possible LECRE origin of the 
6.4 keV line and hard X-ray continuum emission of Sgr 
B2, we implemented a model of low energy cosmic ray 
electron interactions in dense matter to predict the X- 
ray spectrum of their emission and compare it to the Sgr 
B2 measured spectrum. We used a steady-state, thick 
target model, in which accelerated electrons with given 
energy spectra are injected at a constant rate into the 
interaction region and produce atomic reactions as they 
slow down to e nergies below the thre sholds of the var- 
ious reactions (jTatischeff et al.l 119981 ). It includes the 
non-thermal X-ray emission from the secondary, knock- 
on electrons, which mainly result from ionization of am- 
bient H2 molecules and He atoms. The required differ- 
ential ionization cross sections are calculate d from the 
relativistic binary encounter dipole theory (jKim et al.l 
I2OGO) . The continuum emission is due to bremsstrahlung 
of both primary and secondary electrons. Calculations 
of the bremsstrahlung spectra are based on lStrong et al.l 
(|200C1() . The line emission results from the filling of inner- 
shell vacancies produced by both primary and secondary 
electrons in ambient atoms. We considered the Ka lines 
of C, N, O, Ne and Mg (although they are strongly ab- 
sorbed at the GC) and the Ka and K/3 lines from Si, S, 
Ar, Ca, Fe and Ni. The electron-induced K-shell ioniza- 
tion c ross sections a re from Quarles semi-empirical for- 
mula (|Quarleslll976l ) 

The energy spectrum of the fast electrons injected in 
the X-ray production region is a power-law, dQ/dt{E) oc 
E^^ between Emin and Emax- The minimum energy Emin 
is set to 1 keV. The spectral index s and the maxi- 
mum energy Emax are free parameters. So is the ambient 
medium abundance. The model has therefore 3 free pa- 
rameters. 

We note that photoelectric absorption of 
bremsstrahlung X-rays in the source also con- 
tributes to fluorescence line emission. The equiv- 
alent width of the fluorescence Fe Ka line can 
be estimated from the absorbing hydrogen col- 
umn density of the source, Nh (source), to be 
EW = 0.07 X abund x (NH(source)/1023cm-2) keV. 

This approximation is vah d for Nh (source) < lO^'^cm"^ 
(|Leahv fc CreightonI 119931 ). 

The model has been implemented as an XSPEC table 
model and tested on our broadband spectrum of Sgr B2. 
We used the following model: wabs x (apec -I- wabs x 
lecre). To limit the number of free parameters, we 
fixed the temperature of the plasma a t 6.5 keV and 
its metallicity to sola r abundances (see iKovama et al.l 
I2007t iRvuet al.ll2009l ). We obtain a satisfactory fit to 
the XMM-Newton and ISGRI data with a reduced x^ 
of 1.16 (1383 d.o.f.). The CR electron spectral index 
is not well constrained but it is found to be hard, with 
s ^ 1.5, compared to typical acceleration processes. The 
metallicit y of t he molecular cloud relative to solar (as 
given by iLodders, ,2003 ) is A = 3.1 ± 0.2 which is sig- 



nificantly less than that estimated by iRevnivtsev et aLl 
(|2004D . but still too large compared to current measure- 
ments. Thanks to the spectral coverage, the fit provides a 
constraint on the maximum energy of the electrons such 
that Emax = 200 ± 50 keV. This shows that a single pop- 
ulation cannot be at the origin of t he hard X-ray and TeV 
emissions, as has already been bv iCrocker et al.l (|2007D . 
The CR luminosity required to sustain the observed flux 
is 8.9 ± 1.3 X 10'^^ erg/s. This luminos ity requirement is 
slightly lower than that estimated by IRevnivtsev et all 
(J2004) but is still ^ 30% of the total far infrared lu- 
minosity of Sgr B2 (|Goicoechea et afl 120041 ) . We stress 
that since most of this power is dissipated by electrons 
of energy close to E^ax the hard X-ray spectrum is 
important to constrain this parameter. The molecular 
cloud absorption column density is found to be 5.8 ± 
0.3 X 10^'^ cm~^. Th e maximal column density estimated 
bv lProtheroe et aP (.2008) is 12 xlO'^'^ cm~^, suggesting 
that the cosmic ray production and interaction would 
h ave to take place de ep inside Sgr B2. 

iGiveon et al.l ()2002D have measured the GC interstellar 
medium metallicity to be twice solar. Similarly, studies 
of several red giant and supergiant stars in the GC have 
shown that their rnetallicity is only sl ightly above solar 
(jCunha et al.ll2008t [Dalies et al.ll2009D . The metallicity 
required to reproduce the iron line equivalent width is 
therefore too large. 

Besides, the huge power input is a strong limitation, 
since it is close to the total bolometric luminosity of the 
cloud. The decay time is also a strong limitation of this 
scenario: the size of the 6.4 keV line emission region as- 
sociated with the hard X-ray continuum we measure im- 
plies that either the subrelativisitic electron production 
is distributed inside the cloud or that a huge number of 
them have been produced in a limited region and diffuse 
into the cloud. Both solutions are ruled out: the first one 
because we should not expect correlated variation of the 
electron sources all over the cloud, and the second one 
because the time required for electrons to diffuse in a re- 
gion as large as that of the diffuse 6.4 keV line emission 
is much larger than the measured decay time of the sig- 
nal. We can therefore exclude LECRE emission as the 
dominant contribution of the iron line and hard X-ray 
emission in Sgr B2. 

4.2. Hard X-ray reflection from an external distant 

source 

Given the difficulties outlined above with cosmic-ray 
models, it is necessary to turn to the alternative in- 
terpretation in terms of X-ray reflection, for which 
variability on relatively short timescales is expected 
(|Sunvaev fc Ch urazov 1998). When the illuminating 
source is switched off, the decay of the scattered flux 
should occur on a timescale of the light crossing time for 
the cloud. The actual location of the illuminating source 
has then to be determined. An internal source might 
also create such a decrea se if it stopped einitting in the 
recent past. Nevertheless. IRevnivtsev et al.l (|2004l ) reject 
this possibility due to the relative stability of the mea- 
sured Fe Ka flux between 1993 and 2000, which would 
imply that the source was still active in the last decade. 
Since the column density can reach 10^^ cm"'^, we might 
expect an opacity of about 0.66 - not large enough to 
prevent the direct detection of the source during the var- 



Terrier et al. 



ious observations of the molecular cloud in this period. 
It is therefore likely that the source is external. 

We developed an XRN model assuming a relatively 
si mple geometry for the cloud. Following the approach 
of I Mur akami et al.' (2000), we use a 45 pc diameter cylin- 
drical cloud with a symmetry axis parallel to the line of 
sight. The primary irradiating X-ray source is supposed 
to be located at the side of the cloud, at an angle of 90° 
with respect to the line of sight. The region enclosing 
the cloud is divided for the numerical simulations into 
225x225x113 cubic cells of 0.2 pc on a side. 

The density d i stribu tion is assumed to be that of 
iLis &: GoldsmithI ()1990D . with a very high density 5-pc- 
wide core surrounded by a more diffuse envelope. The 
total mass of the cloud is normalized to be 1.9 10^ M0, 
w hich is the virial mass o f Sgr B2 as recently estimated 
bv lProtheroe et all ('2008'). 

The probability of scattering of a primary X-ray pho- 
ton is calculated from both the coherent (Rayleigh) and 
incoherent (Compton) scattering cross sections, taking 
into account the 13 most abundant elements in the am- 
bient medium (H, He, C, N, O, Ne, Mg, Si, S, Ar, 
Ca, Fe and Ni). These cross sections ar e taken from 
the X COM photon cross sections database (jBerger et al.l 
I2009D . Above a few keV depending on the composition, 
the total scattering cross section is dominated by the in- 
coherent process. The energy of the scattered photon is 
then given by E'^ = ExrrieC^ / {Ex + rrieC^) {Ex being 
the energy of the primary photon and ?Tie the electron 
mass), which is appropriate for a Compton scattering at 
90°. 

The fluorescent Ka and K/3 lines of Si, S, Ar, Ca, Cr, 
Mn, Fe and Ni are considered. We also include the Ka 
line of C, N, O, Ne and Mg although these lines are 
strongly absorbed in the present case. The K-shell photo- 
ioniza tion cross sections are from IVerner fc YakovlevI 
(|1995D . Photoelectric absorption of both the primary and 
secondary X-rays is calculated from the elem e ntal c ross 
sections of iBalucinska-Church fc McCammonI ()1992f) . It 
thus depends on the cloud composition. 

The resulting spectra were then implemented in an 
XSPEC table model with 3 parameters: the irradiating 
source spectral index (which can vary from 1 to 5), the 
ambient medium abundance, which can vary fro m 0.5 to 
3 time s solar (the solar metallicity is taken from lLoddersI 
(|2003| )) and the flux normalization. 

We used a total model similar to that used in the 
previous section: wabs x (apec -I- xrn). We obtain a 
satisfactory fit to the XMM-Newton (2004) and ISGRI 
data (2003-2004) with a reduced x^ of 1.16 (1384 d.o.f.). 
The illuminating source power law spectral index is 
found to be r ^ 2, the abundance relative to solar 
A = 1.3 ± 0.1 and the average luminosity of the illu- 
minating source required to sustain the observed flux is 
L(l - 100 keV) = 1.1 x 10^9 (d/lOOpc)^ erg/s, where d is 
the distance of the source to Sgr B2. 

As has been noted in previous work, if the illuminating 
source were to be at a distance of ~ 10 pc its luminosity 
would have been larger than 10^^ erg/s for several years 
and it should have been detected in the past decades, but 
if the source is farther away the luminosity constraints 
leave only Sgr A* as a viable candid ate ( Kovama ct al. 
[l996t iMurakami et al.l[2000l : iRevnivtsev et al.ll2004D . In 



the latter case, Sgr A* should also be responsible for 
the measured Fe Ka in some other molecular clouds of 
the CMZ. The observations of time dependent effects in 
other places is therefore very strong evidence of the XRN 
nature of the Sgr B2 hard X-ray emission a nd supports 
Sgr A * as the origin of the illuminating flux (jPonti et al.l 
I20T0I) . 

In a recent analysis of VLBI observations, Reid et al. 
(2009) have measured the parallax of the Sgr B2 (M) 
and Sgr B2 (N) dust condensates. The authors have es- 
timated the Galactocentric radius at Rq = 7.91q7 kpc 
and the proper motion of Sgr B2 relative to Sgr A* sug- 
gests that the former is 130 ± 60 pc closer to us than 
the latter. This fact is also suggested by the study of 
photoelectric absorption in front of th e variou s molecu- 
lar clouds inside the Sgr B complex ()Rvu et a l. 2009). 
In a scenario where intense Sgr A* activity is the origin 
of the hard X-rays from the cloud one can date the end 
of the high state period to be IOOI25 years ago, notably 
different from the widely discussed 300 years obtained 
assumin g the GMC is at the same distance as Sgr A* 
(see e.g. IKovama et al.l (|1996( )). 

This geometry is inconsistent with the model we use 
which requires orthogonal photon incidence on the cloud. 
At 130 pc in front of Sgr A*, several effects should modify 
the resulting spectrum. First, the mean column density 
seen by scattered photons is larger, increasing the absorp- 
tion with respect to the orthogonal scenario. Second, the 
Compton diffusion cross section is larger for scattering 
angles smaller than 90° at high energies. The resulting 
spectrum should therefore be harder and brighter in the 
hard X-ray domain and more attenuated under < 15 keV. 
To estimate the effect on the measured spectral parame- 
ters, we have performed simple Monte-Carlo simulations 
of photon propagation in a spherical uniform cloud. For 
clouds mass comparable to that of Sgr B2, we found out 
that the differences on the measured spectral index could 
be up to 0.2 in the 15-100 keV range for cos 6* ranging 
from to 0.9 (corresponding to a cloud 190 pc in front of 
the GC). The total scattered flux also increases by up to 
30% in the 5 to 100 keV band. Therefore, taking into ac- 
count uncertainties due to the relative position of SgrB2 
and Sgr A*, the actual luminosity from the flare should 
range between 1.5 and 5 10^^ erg/s and its spectral index 
should be r = 2 ± 0.2. 

5. CONCLUSIONS 

INTEGRAL has been regularly monitoring the Galac- 
tic Centre for 7 years. U sing the complete dataset, we 
confirm the result from IRevnivtsev et al] ()2004D that 
IGR J17475-2822 corresponds to hard X-ray (>20 keV) 
emission from the Sgr B2 molecular cloud. In partic- 
ular, it is in good positional agreement with the very 
massive and dense core of the cloud. The larger and 
more diffuse component might be contributing but its 
flux should be diluted because of its extent relative to 
the PSF size. Since the mean density is of the order of 
1 — 3 X lO'^cm^^, the envelope column density is smaller 
by a factor of a few compared to the core column density 
(which is close to lO^'^cm^^), so that the diffuse hard X- 
ray contribution has to be much smaller than that of the 
core. 

The regular monitoring shows that the Sgr B2 20-60 
keV emission is fading, with a flux decrease of about 40% 
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Fig. 5. — Broadband spectrum of the Sgr B2 inoleeular cloud fitted with a X-ray nebula spectral model (see text for details) 



after 7 years of observations. A test of hypothesis shows 
that a hnear decay is favored compared to a constant 
flux at the 4.8 a confidence leveL The hnear decay char- 
acteristic timescale is found to be ~ 8 yrs, comparable 
with the light-travel time across the cloud core. This 
result is obtained thanks to the very deep exposure ob- 
tained by INTEGRAL on the GC regions and does not 
suffer from instrument to instrument intercalibration er- 
rors. It is consistent with the results obtained comparing 
the 6.4 keV fluxes measured by Chandra, XMM -Newton 
and Suzaku from 2000 to 2005 (llnui et al.|[2009l) . 

Cosmic ray models of the iron line and continuum emis- 
sion are put to a test by this fast variability. Inverse 
bremsstrahlung emission from subrelativistic ions is ruled 
out because of the long cooling timescale of these parti- 
cles in the dense environment of the cloud core. Low 
energy cosmic ray electrons are not formally ruled out: 
the lifetime of '^ 100 keV electrons in the dense matter 
of Sgr B2 is smaller than the measured decay time, so 
that if the cosmic-ray injection had ceased with a simi- 
lar or smaller timescale, we might obtain a similar flux 
decrease. Nevertheless, combining the measured decay 
with the large electron injection power required (which 
we fitted to be of the order of 10*" erg/s) we can exclude 
this scenario as well. 

Sgr B2 is therefore most likely a reflect i on ne bula, 
as was initiall y prop osed by iSunvaev et al.l ()1993[ ) and 
iKovama et al.l ()1996D . The location of the illuminating 
source is yet to be determined unequivocally. As was 
already noted by several works, an internal or a nearby 
illuminating source should have been detected in the last 
decade or so. The most likely explanation is therefore 
a period of high state of Sgr A*. The observation of 
variability of the Fe Ka line in other massive molecu- 



lar clouds of the C MZ strongly supports this hypothesis 
(jPonti et al.l 120101 ) ■ since several distinct bright illumi- 
nating sources active for at least several years would be 
required to light all these clouds. They could be illumi- 
nated by the same outburst from Sgr A*. 

Usi ng a m o del based on the approach of 
Mura kami et al.l (|2000[ ). we tested this scenario and 
found that the illuminating spectrum must be hard 
(photon index ~ 2), which is well in the range of mea- 
sured Low Luminos ity Active Ga.lactic Nuclei (LLAGN) 
spect ral indices (|Gu fc Cad l2009t iPellegrini et al.l 
l2000f ). Since parallax measuren ients place Sgr B 2 130 
± 60 pc in front of Sgr A* (jReid et all l2009l ). our 
measurement implies a mean luminosity of 2-5 xlO^^ 
erg/s, an increase of more than 5 orders of magnitude 
compared to the current one placin g it among faint 
LLAGNs (IHoI [2001 IGu fc Cad l200l . At the end of 
the outburst the variability of Sgr A* is constrained 
to be faster than the measured decay rate, which is 
compar able to the long term variab il ity observed in 
M81* (|Ivomoto fc Makishima] 1200X1 : iLa Parola et al.l 
120041 ). The end of the episode of high luminosity in 
Sgr A*, which is currently being mirrored by Sgr B2, 
occurred about 100 yrs ago, with the intervening period 
seeing the output of the central source fall by some 
5 orders of magnitude. Assum ing that Sgr B2 is at 
the same distance that Sgr A*, llnui et ahl (|2009[ ) have 
proposed that the luminosity of the latter has been 
continuously decreasing with a constant decay time 
for 300 years. The new position implies that Sgr A"^ 
luminosity had to decrease on a much shorter timescale. 
The origin of such large swings in the luminosity of the 
supermassive black hole at the centre of our Galaxy is 
presently unknown but the subject of intensive on-going 
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